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ABSTRACT 

We use the Monte-Carlo Markov Chain method to explore the dark energy 
property and the cosmic curvature by fitting two popular dark energy parame- 
terizations to the observational data. The new 182 gold supernova la data and 
the ESSENCE data both give good constraint on the DE parameters and the 
cosmic curvature for the dark energy model wq + w a z/(l + z). The cosmic cur- 
vature is found to be \ flk\ < 0.03. For the dark energy model wo + w a z/ (1 + z) 2 , 
the ESSENCE data gives better constraint on the cosmic curvature and we get 

\n k \ < 0.02. 



Subject headings: Cosmology: cosmological parameters — Cosmology: observa- 
tions 
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Introduction 



The supernova (SN) la observations indicate the accelerated expansion of the Universe 



( IRiess et al.lll998t iPerlmutter et al.lll999l ). The direct and model independ ent evidence of 



the acc elerat ion of the Universe was shown by using the energy conditions in IGong k Wang 



( l2007bl ) and IGong et al.l (120071 ). The driving force of the late time acceleration of the Uni- 
verse, dubbed "dark energy (DE)", imposes a big challenge to theoretical physics. Although 
the cosmological constant is the simplest candidate of DE and consistent with current ob- 
servations, other possibilities are also explored due to many orders of magnitude discrep- 



ical constant. For a review of DE models, s e e for exampl e 
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Due to the degeneracies among the parameters in the model, complementary cosmological 
observations are needed to break the degeneracies. The Wilkinson Microwave Anisotropic 
Probe (WMAP) measurement on the Cosmic Microwave Background (CMB) anisotropy, to- 
gether with the SN la observations provide complement ary data. In this paper, we use the 



three-year WMAP (WMAP3) data 



Wood-Vasev et al.ll2007l : iDavis et al. 



Spergel et all 120071) . the SN la data (IRiess et al.l 12006 



20071 ) a nd the Baryon Acousti c Oscillation (BAO) mea- 



surement from the Sloan Digital Sky Survey (jEisenstein et al.ll2005l ) to study the property of 
DE and the cosmic curvature. Two DE models w(z) = Wn+w n z/(l+z) dChevallier k D. Polarski 



2001 



Linder 2003) and w(z ) = wq + w a z/(l + z) 2 ( Jassal et al. 20051 ) are considered. In 



Elgar0y k Multamaki (j2007 ) . the authors showed that combining the shift parameters R 
and the angular scale l a of the sound horizon at recombination appears to be a good approx- 
imation of the full WMAP3 data. Wang and Mukherjee gave model independent constraints 
on R and l a by using the WMAP3 data, they also provided the covariance matrix of the 
parameters R, l a and VL^h 2 (IWang k Mukherj eel 120071 ). So we use the shift parameter R, the 
angular scale l„. of the sou nd horizon at recombination and their covariance matrix given in 
Wang k Mukherj ed (120071 ) instead to avoid using several inflationary model parameters and 
calculating the power spectrum. When the covariance matrix is used, we have six parameters. 
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We use the Monte-Carlo Markov Chain (MCMC) method to explore the p arameter space. 



Our MCMC code is based on the publicly available package COSMOMC (ILewis &: Bridle 



20021 ). 



The paper is organized as follows. In section II, we give all the formulae and show the 
constraint on the cosmic curvature is much better by using the parameters R, l a and their 
covariance matrix than that by using the parameter R only. We also discuss the effect of the 
radiation component fl r on l a . In section III, we give our results. We discuss the analytical 
marginalization over H Q in appendix A. 



2. Method 



For the SN la data, we calculate 



X 



a? 



where the extinction-corrected distance modulus n(z) = 51og 10 [dz,(z)/Mpc] + 25, Oi is the 
total uncertainty in the SN la data, and the luminosity distance is 

dz' 



dr.(z) 



l + z 



:S1I1I1 



E z' 



(2) 



here 



sinn 



sm(y/\Q k \x)/ y/\Q k 



x. 



if VL k < 0, 
if fife = 0, 

if n k > o, 



and the dimensionless Hubble parameter is 

E 2 {z) = H 2 (z)/Hl = Q m (l + zf + Q r (l + zf + Q k (l + zf + Q DE , 



(3) 



(4) 



where Q = SirGp/ (3Hq), p r = a^T^ mh , is the Stefan-Boltzmann constant, the CMB 
temperature T cm b = 2.726K, and Qde is the DE density. Note that the distance normal- 
ization is arbitrary in the SN la data, the Hubble constant H determined from the SN 
data is also an arbitrary number, not the observed Hubble constant. Therefore we need to 
marginalize over this nuisance parameter Hq. The parameter Hq is marginalized over with 
flat pri or, the analytical marKinalizat i on method i s discussed in Appendix A. For the DE 
model Jchevallier fc D. PolarsklboOll ; lLinderil2003h 

w a z 



w(z) = w + 



l + z' 



(5) 
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the dimensionless DE density is 

Vde(z) = (1 - n m - n k - n r ){l + z ) 3( - 1+Wo+w ^ exp[-3 Wa z/(l + z)]. 



For the DE model (j.Tassal et all 120051 ) 



w(z) =w + 

(l + zy 

the dimensionless DE density is 

n DE {z) = (1 - n m - V k - fi r )(l + z)^ 1+wo) exp [3™ z 2 /2(l + zf] . 

For the SDSS data, we add the term 

"A- 0.469(0.95/0.98)-°- 35 " 
0.017 



(6) 



(7) 



to x 2 ( lEisenstein et al.ll2005l ; ISpergel et al.l 120071 ). where the BAO parameter 
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For WMAP3 data, we first add the term 



R- 1.71 
0.03 



to x 2 ( IWang fc Mukherjed 120071 ). where the shift parameter 



R 



:smn 



dz 
E(z) 



(10) 



and z is = 1089 ± 1. 



When we fit the DE models (jSJ) and (j7]) to the observational data, we have four param- 
eters Q m , Qk, w and w a . The MCMC method is used to explore the parameter space. The 
marginalized probability of Q k is shown in Fig. [TJ It is obvious that the cosmic curvatu r e can- 
not be well constrained f or the DE model (JHJ). As discussed in lElgarOy fc Multamakil (120071 ) 
and I Wang fc Mukherjed (120071 ) , the combination of the shift parameter and the angular scale 
of the sound horizon at recombination gives much better constraints on c osmological parame- 
ters. So we add the angular scale of the sound horizon at recombination (IWang fc Mukherjee 
2007h 

TlR/y/n^ 



!Zdzc s /E{z) 



302.5 ± 1.2, 



(11) 
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where the sound speed c s = 1/ y 3(1 + Rha), Rh = 3l5000£lhh (T rmh /2.7K) , a is the scale 
factor, and Q h h 2 = 0.02173±0.00082 jWang fa Mukherjeell2M7l ). To implement the WMAP3 



data, we need to add three fitting parameters R, l a and Q h 2 . So we need to add the 
term Ax,;Cov~ (xi,Xj) Ax j to x 2 > where xi = (R, l a , Vt^h 2 ) denote the three parameters 

and Cov(xj,Xj) is the covariance matrix for the three 
3, we use the covariance matrix for Xi = (R, l a , Q h 2 ) 



for WMAP3 data, Ax,- 



X. ,; 



x 



obs 



derived in 


Wane fa Mukheriee 


Wang fa Mukherjee 


( 


2007 


) 15 



element in the matrix is obtained by marginalizing over all other variables. Therefore, the 
covarianc e matrix between Xj and xj is the three by three sub-matrix of the full six by six 
matrix in I Wang fa Mukherjed (120071 ). The marginalized probability of Q%. is shown in Fig. 
[U We see that the cosmic curvature is constrained better with the addition of the angular 
scale l a of the sound horizon at recombination. 

Since the angular scale of the sound horizon depends on the early history of the Universe, 
so it strongly depends on Q r . However, we can neglect the effect of Q r when we evaluate the 
distance modules p(z) and the shift parameter R because the Universe is matter dominated. 
So only when we implement the CMB data with l a , we need to consider the effect of Q r . 
We know the energy density p r of radiation, so the dependence of fl r = 87tGp r / (3ifg ) is 
manifested by the Hubble constant H . Since we can neglect the effect of Q r in fitting 
SN la data, so the effect of the observed value of H can be neglected by marginalizing 
over it. Therefore, we use the Hubble constant Hq as a free parameter instead of Q r . The 
marginalized probabilities of Qk for H = 65 km/s/Mpc and H = 72 k m/s/Mpc are shown in 
Fig. [21 We see that the results indeed depend on H . As discussed in ( lElgarpy fa Multamaki 



20071 ). the combination of R and l a approximates the WMAP3 data and the WMAP3 data 
depends on H through l a . So, as expected, l a also depends on H . From now on we also 



take H q as a fitting parameter, and impose a prior of H = 72 ± 8 km/s/Mpc (IFreedman 



200ll ). To understand why we can marginalize over Hq in fitting SN la data and treat Hq as 
a parameter in fitting WMAP3 data, we should think that we actually treat Q r , not Ho as 
a parameter when fitting the WMAP3 data. The parameter H is not the observed Hubble 
constant when fitting the SN data because the normalization of the distance modulus was 
chosen arbitrarily. In summary, we have six fitting parameters for the DE models (jSJ) and 
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Results 



In this section, we present our res ults. We first use the 182 gold SN la data (IRiess et al. 



2006), then we use the ESSENCE data (IRiess et al.ll2006l : lWood-Vasev et al.ll2007l : lDavis et al 
20071 ). F or the SN la data, we consider both the SN la flux averaging; w ith marginalization 
over H ( IWanglboOOl ; IWang fc Mukherj eel 12004 ; IWang fc Tegmarkl 12004 ) and the analytical 
marginalization without the flux averaging. The results with the analytical marginalization 
are shown in solid lines and the results with flux averaging are shown in dashed lines. We 
also put the ACDM model with the symbol + in the contour plot. 



3.1. Gold SN la data 

Fig. [3] shows the marginalized probabilities for Q m , Q k , w and w a for the DE model w + 
w a z/(l + z). Fig. H] shows the marginalized fl m -fik a nd Wo-w n contours. The W o-w a contour 



with the flux averaging is consistent with the result in lWang fc Mukherjed ( 120071 ). From Figs. 
[3]and|H we see that the difference in the results between the analytical marginalization and 
the flux averaging is small. The ACDM model is consistent with the observation at the la 
level. The value of w a is better constrained with the analytical marginalization. 

Fig. [5] shows the marginalized probabilities for Q m , fl k , wo and w a for the DE model 
Wo + w a z/{l + z) 2 . Fig. [6] shows the marginalized Q m -Qk and w -w a contours. From Figs. 
E]and [HI we see that the parameters are a little better constrained with the flux averaging. 
For the analytical marginalization, the ACDM model is consistent with the observation at 
the 2a level. For the flux averaging, the ACDM model is consistent with the observation at 
the la level. 



3.2. ESSENCE data 

Fig. [7] shows the marginalized probabilities for Q m , fl k , w and w a for the DE model 
wo + w a z/{l + z). Fig. [8] shows the marginalized fl m -flk and wo-w a contours. From Figs. [7] 
and [HJ we see that the difference in the results between the analytical marginalization and 
the flux averaging is small. The ACDM model is consistent with the observation at the la 
level. 

Fig. [9] shows the marginalized probabilities for Q m , Qk, wq and w a for the DE model 
wq + w a z/(l + z) 2 . Fig. [TUI shows the marginalized Q m -Qk and wo-w a contours. From Figs. 
M and [TUI we see that the parameters are a little better constrained with the analytical 
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marginalization. The ACDM model is consistent with the observation at the la level. 

We summarize the results in Tables 1 and 2. We do not see much improvement on 
the constraints on the DE parameters and the cosmic curvature by using the flux averaging 
method. For the DE model wq + w a z/(l + z), the gold data gives better constraints than the 
ESSENCE data on the DE parameters wq and w a , but both data give good constraints on 
the cosmic curvature. For the DE model wo + w a z/ (1 + z) 2 , the ESSENCE data gives much 
better constraint on the cosmic curvature than the gold data, although the constraints on the 
DE parameters w and w a are almost the same for both data. For the 182 gold data, the DE 
model wo + w a z/(l + z) gives much better constraints on the cosmic curvature flk- For the 
ESSENCE data, the two DE models give almost the same constraint on Q m and For the 
DE model Wq + w a z/ (1 + z), the mean value of Wq determined from the observation tends to 
be Wq > —1, while the mean value of Wo is less than —1 for the DE model Wq + w a zj (1 + z) 2 . 

From Tables 1 and 2, we see that the constraints on Qk are almost the same for the 
two different DE models ((Sj) and ((Tj). In other words, the resu lts we obtained on flk do not 



depend on the chosen models much. Recently, the authors in IClarkson et al.l (120071 ) found 
that the assumption of a flat universe induces critically large errors in reconstructing the 
dark energy equation of state at z > 0.9 even if the true cosmic curvature is very small, 
Qk ~ 0.01 or less. They obtained the result by fitting the data derived from a DE model 
with Qk 7^ with a flat model, so the result may not be conclusive. To see how the value of 
Qk affect the constraints on the property of DE, we perform the MCMC analysis on the DE 
models (jSJ) and (j7|) with flk = 0. The results are reported in Tables 3 and 4. Although the 
uncertainties of Qk change the values of wo and w a , the ranges of wo and w a are almost the 
same for small flk- 

In conclusion, we first confirm previous results that the shift parameter R alone does 
not give good constraint on Q k , we must combine R and l a to constrain Q k - By using R, l a 
and their covariance matrix, we get almost the same results as those obtained by using the 
original WMAP3 data. Without calculating the power spectrum, the fitting process is much 
faster and efficient. The cosmic curvature is found to be |QJ < 0.03. 



YGG and AW thank Yun Wang for the help with the MCMC method. YGG is grateful 
of Zong-hong Zhu for fruitful discussions, and he is supported by NNSFC under grant No. 
10605042. A. Wang is partially supported by the VPR funds, Baylor University. 
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A. Analytical marginalization on Hq 

By assuming a flat prior P(Hq) = 1 for Hq, the marginalization over Hq means 

L = e~ x ™ /2 = J e~ x2/2 P(H )dH = J e - x2/2 dH . (Al) 

Let x = 51og 10 if and a { = ii obs {zi) - 25 - 5 log 10 [(l + Zi)sian(^/\Q k \ dz' /E(z'))/ ^\Q k \], 
and substitute Eq. ([T]) into the above Eq. flAll) . we get 



L = ^2/ dxexp 



^ Jcixexp 



1 (Qi+a:) z InlO™ 

I" — X 



-i(e^) ( 

llEj ai/^-lnlO/5) 2 



In 10 

5 



1/2 



l(Ei«i/g ?-lnlO/5) a - 

W 1 ,. ■ 



So the minimum \ is 



(A2) 



y. o? _ (E,^/j-lnl0/5) 2 _ 2 ^ / In 10 



2tt 



(A3) 
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Table 1: The marginalized results with la errors for the model Wq + w a z/(l + z) 
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Analytical 
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-0.94 ±0.25 




W a 


0.34 ±0.77 


-0 0S +L04 


-0.70±l;lf 


-0.26ii;|? 




Fig. 1. — The marginalized probabilities of f^. The solid lines denote the results using the 
shift parameter R, the angular scale l a , and the full covariance matrix. The dashed lines 
denote the results using the shift parameter only. The black lines are for the dark energy 
model Wq ± w a z/(l ± z) 2 and the red lines are for the model wq ± w a z/ (1 ± z). 



-0.02 0.02 0.04 0.06 0.08 0.1 0.12 0.14 

Fig. 2. — The marginalized probabilities of Q]~. The solid lines denote the results with 
Hq = 65. The dashed lines denote the results with Hq = 72. The black lines are for the dark 
energy model wq + w a z/ (1 + z) 2 and the red lines are for the model u> + w a z/ (1 + z). 



Table 2: The marginalized results with lcr errors for the model Wq + w a z/(l + z) 2 
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Fig. 3. — The marginalized probabilities for the DE model Wq + w a z/(l + z) by using the 
gold SN la data. The solid lines denote the results with analytical marginalization and the 
dashed lines denote the results with flux averaging. 



Table 3: The marginalized results with la errors for the model Wo + w a z/(l +z) with Qk = 
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Fig. 4. — The marginalized la and 2a Q m -Q k and w -w a contours for the DE model w + 
w a z/(l+z) by using the gold SN la data. The upper panels denote the results with analytical 
marginalization and the lower panels denote the results with flux averaging. 



Table 4: The marginalized results with la errors for the model w + w a z/(l+z) 2 with Q k = 
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Fig. 5. — The marginalized probabilities for the DE model w + w a z/(l + z) 2 by using the 
gold SN la data. The solid lines denote the results with analytical marginalization and the 
dashed lines denote the results with flux averaging. 



-16- 




Fig. 6. — The marginalized la and 2a Q m -Qk and WQ-w a contours for the DE model Wq + 
w a z/(l + z) 2 by using the gold SN la data. The upper panels denote the results with 
analytical marginalization and the lower panels denote the results with flux averaging. 
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Fig. 7. — The marginalized probability distributions for the dark energy model wo+w a z/ (1 + 
z) by using the ESSENCE data. The solid lines denote the results without flux average and 
the dashed lines denote the results with flux average. 
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Fig. 8. — The marginalized Vt m -Vtk and WQ-w a contours for the dark energy model Wq + 
w a z/(l + z) by using the ESSENCE data. The upper panels denote the results without flux 
average and the lower panels denote the results with flux average. 
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Fig. 9. — The marginalized probability distributions for the dark energy model wo+w a z/ (1 + 
z) 2 by using the ESSENCE data. The solid lines denote the results without flux average and 
the dashed lines denote the results with flux average. 
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Fig. 10. — The marginalized Vt m -Vtk and wo-w a contours for the dark energy model wq + 
w a z/ (1 + z) 2 by using the ESSENCE data. The upper panels denote the results without flux 
average and the lower panels denote the results with flux average. 



